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ABSTRACT 

The effects of rapid rotation on stellar evolution can be profound. We are now 
beginning to gather enough data to allow a realistic comparison between different 
physical models. Two key tests for any theory of stellar rotation are first whether it 
can match observations of the enrichment of nitrogen, and potentially other elements, 
in clusters containing rapid rotators and secondly whether it can reproduce the ob- 
served broadening of the main sequence in the Hertzsprung-Russel diagram. Models 
of stellar rotation have been steadily increasing in number and complexity over the 
past two decades but the lack of data makes it difficult to determine whether such 
additions actually give a closer reflection of reality. One of the most poorly explored 
features of stellar rotation models is the treatment of angular momentum transport 
within convective zones. If we treat the core as having uniform specific angular mo- 
mentum the angular momentum distribution in the star, for a given surface rotation, is 
dramatically different from what it is when we assume the star rotates as a solid body. 
The uniform specific angular momentum also generates strong shears which can drive 
additional transport of chemical elements close to the boundary of a convection zone. 
A comparison of different models and their reproduction of observable properties with 
otherwise identical input physics is essential to properly distinguish between them. We 
compare detailed grids of stellar evolution tracks of intermediate and high-mass stars 
produced using several models for rotation generated with our new stellar rotation 
code. 

Key words: stars:evolution, stars:general, stars:rotation 



1 INTRODUCTION 

It has been known for many years that rapid rotation can 
cause significant changes in the evolution of stars (e.g. 
Kippcnha hn. Mever-Hofmeister fe Thom as 1970). Not only 
does it cause a broadening of the main sequence in the HR 
diagram but it also produces enrichment of a number of 
different elements at the stella r surface (e.g. iHunter et al.l 
2009.; iFrischknecht" et al. 2010). With new large scale sur- 
veys, such as the VLT-FLAMES survey of massive stars, 
now reaching maturity, t he data available for rota ting stars 
are growing rapidly (e.g. lEvans et al.l 12005. 2006). Any vi- 
able model of stellar rotation must be able to match the ob- 
served changes in chemical enrichment and structure. The 
treatment of rotation and its induced chemical mixing in 
stars has cha nged dramat ically over the past two decades. 
The model o f lZahrJ (| 19921 ) has formed the basis for much of 
this work and many variations from the original have been 
used during this tim e to generate different sets of stellar 
evolution tracks (e.g. iTalon et "all 1 19971 ; iMevnet fe MaedeJ 
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l200d; lMaederll2003l). Altern a tive f ormalisms, such as that 



of Heger. Langer fc Wooslevl |2000l ). treat the physical pro- 
cesses very differently. Particular emphasis is often placed on 
the treatment of meridional circulation. While those models 
based on Zahn's (1992) treat meridional circulatio n as an 
advective process, iHeger. Langer fc Wooslevl (gOOO) treat it 
as a diffusive process. This leads to fundamentally different 
behaviour. This is just one feature of the model which may 
lead to significantly different results. Even so, both treat- 
ments are used and frequently quoted in the literature for 
their predictions of the effect of rotation on stellar evolution. 

One of the most poorly explored features of stellar rota- 
tion is the treatment of angular momentum transport within 
convective zones. Current ID models, treat convective zones 
as rotating solid bodies. This is not necessarily correct and 
there are strong reasons to explore alt ernatives such as 
unifo r m specific angular m omentum (e.g. lArnett fc Maekir] 
120091 ; iLesaffre et alj |2010|) . This potentially has dramatic 
consequences for the evolution because a star with uniform 
specific angular momentum through its convective core has 
much more angular momentum for a given surface rotation 
velocity. It also has a strong shear layer at the convective 
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boundary that can drive additional transport of chemical 
elements. 

Different models for stellar rotation have not been com- 
pared directly on a common numerical platform alongside 
otherwise identical input physi cs before. From the Cam- 
bridg e stellar evolution code |Eggletonl Il97ll ; I Pols et all 
Il995l ) we have produced a code capable of modelling ro- 
tating stars in ID under the shellular rotation hypothesis of 
IZahni (|l992l ). The code, RoSE (Rotating Stellar Evolution), 
can be easily programmed to run with different physics for 
stellar rotation and can model both radiative and convective 
zones under a range of different assumptions. This allows us 
to compare a variety of models for stellar rotation and de- 
termine what, if any, observable traits could be used to dis- 
tinguish between them. We foresee two possibilities; either 
we can identify clear observational tests to eliminate certain 
models or the models show no testable difference in which 
case a simplified model could be formulated to provide the 
same results. 

In section 2 we outline the physical ingredients to RoSE 
and the different models already implemented. In section 3 
we present a comparison of the evolutionary predictions for 
each model. In section 4 we present our summary and con- 
clusions. 



2 STELLAR EVOLUTION CODE 

Here we present the physical ingredients of the numerical 
code, RoSE, used to implement the models of stellar rota- 
tion. 



2.1 Stellar evolution package 

RoSE is based on the Cam bridge ste l lar ev olution code, 
STARS. Originally written bv lEggletor] (|l97ll ), the code has 
been modified many times over the pa st thirty years. For de- 
tails o f the last significant update see IStancliffe fc Eldridge! 
(2009). We now solve the four structure equations, seven 
chemical equations and now the angular velocity equation in 
a single, implicit, Newton-Raphson iterative step. We calcu- 
late 1 H, 4 He, 3 He, 12 C, 14 N, 16 and 20 Ne implicitly and 39 
other isotopic abundances can be calcula ted explicit l y. The 
nucl ear reaction rat e s wer e updated by iPols et all (|lQ95l l 
and IStancliffe et al. I (I2005T) and are based on the reaction 
rates of lCaughlan fc Fowled (Il988l'l. The opacities were last 
updated by (Eldridge fc Tout! J2004h and are based on the 
OPAL op acities (llglesias fc Rogers! 19961 ') at high tempera- 
tures and iFerguson et all (|2005l ) at low temperatures. Also 
inclu ded are further corrections for major molecular opaci- 
ties (Stancliffc & Glcbbcck 2008). The equation of state is 
that of lPols et all (Il995l). convectio n is treated by mixing- 
length theory ( Bohni- Vitensd 1958TI and a model fo r con- 
vective overshooting ( Schroder. Pols fc Egglcton 1997) is in- 
cluded. Whilst the code has a number of models for mass 
lo ss programmed, we restrict oursel ves to the mass-loss rates 
of lVink. de Koter fc Lamersl (120011 ) for massive stars. These 
rates apply to non-rotating stars and are modified as ex- 
plained in section [2T4l 



2.2 Structure equations for rotating stars 

The centrifugal force caused by rotation affects the hydro- 
static balance of the star, effectively reducing the local grav- 
ity. On a surface of constant radius the centrifugal force acts 
more strongly at the equator than the poles so the distor- 
tion of the star depends on co-latitud e and our assum ption of 
spherical symmetry is no longer valid. lTassoull|l978h showed 
that, except for stars that are close to critical rotation, the 
effect of rotational deformation remains axially symmetric. 
Enhanced mass loss from the surface because of rotation gen- 
erally keeps stars rotating sufficiently below critical. Even 
when this is not the case it is only the outer most layers that 
are affected. In models where the angular momentum dis- 
tribution in convective regions is uniform the rotation rate 
may approach critical there but because convective turbu- 
lence is already considered to be fully asymmetric, we do 
not need to consider further non-axial instabilities owing to 
the rotation. 

We adopt similar adjustme nts to the stellar stru cture 
equations to those descr ibed by lEndal fc Sofia! (Il976l ) and 
iMevnet fc Maederj (|l997T ). First we define Sp to be a surface 
of constant pressure, P. Vp is the volume contained within 
Sp and rp is the radius of a sphere with volume Vp = 
47rrf>/3. The equation of continuity is then preserved in its 
non rotating form, 



drp 



4irrpp, 



(1) 



where mp is the mass enclosed within Sp and p is the density 
on the isobar which is assumed to be uniform. As we discuss 
in section [231 we expect variables and chemical abundances 
to be uniform across isobars owing to the strong horizon- 
tal turbulence c aused by t he strong density stratification 
present in stars (|Zahnlll992l ). The local gravity vector is 



Soft 



-— + Q, 2 Rsm 2 e ) e r + (fi 2 i?sin(9cos6>) eg, (2) 



where Q is the local angular velocity. To proceed further we 
define the average of a quantity over Sp as 



< q 



(3) 



where da is a surface element of Sp. Using this notation the 
equation of hydrostatic equilibrium becomes 

dP Gm P 



dmp 



471-rf, 



fp, 



where 



fp = 



4irri 



< 9 J > 1 



(4) 



(5) 



GmpSp 

and g c ff = \g cH \- Hence with the new definition of variables 
we can retain the same ID hydrostatic equilibrium equa- 
tion modified by a factor of fp which tends to unity for no 
rotation. 

The equation for radiative equilibrium is similarly mod- 
ified to 

dlnT 3kPL p f T 



d\nP WnacGmpT* f P ' 



(6) 



where Lp is the total energy flux through Sp, P is the pres- 
sure, T is the temperature, k is the opacity, a is the radiation 
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constant, c is the speed of light, G is the gravitational con- 
stant and 



h = VST 



1 ^ N- 1 



(< ffeff >< Seff >) 



(7) 



Again, the non-rotating equation for stellar evolution has 
been preserved except for the multiplication by fr/fp- Of 
the two factors, fp deviates further from unity for a given ro- 
tation than fx- Additional secondary effects of the reduced 
gravity must be taken into account when calculating quan- 
tities such as the pressure scale height and Brunt- Vaisala, 
frequency. For the remainder of the paper we drop the sub- 
script P. 



2.3 Meridional circulation 

The amount of thermal flux, F through a point in a star 
behaves as F oc g e g(9). As in equation ©, g cB is strongly 
dependent on co-latitude and so the radiative flux is greater 
at the poles than at the equator. This leads to a global 
thermal imbalance that drives a meridional circulation. The 
presence of meridiona l circulation has b een considered for 
nearly a century since [Von Zeipell (|l924l ). In the past it has 
com monly been approximated by Eddington-Sweet circula- 
tion |Sweetlll950l ) . IZahnl (| 19921 1 proposed an alternative but 
similar treatment of the meridional circulation based on en- 
ergy conservation along isobars and this is the formulation 
we are currently using. In spherical polar coordinates the 
circulation takes the form 

U = E/(r)P 2 (co S 0)e r + V(r) m ^° S ^ e e , (8) 

a u 

where U and V are linked by the continuity equation 

v = ^4-,( P r 2 u). 



6pr dr 



(9) 



and P2 is the second Legendre polynomial P2 (x) = \ (3a; 2 — 
1) . We current ly use an approximate f orm for the meridional 
circulation bv lMaeder fc Zahnl (|l998f l 

T P 1 

U = Co 



meffgeff CppT Vad - V + V M 

2 \ /„r>2„3\ ( 10 ) 



E m 2nGp ) \ 3Gm 



where m c s = m 



\ L 2ttGp)> 



£ = E n 



+ E K 



the total 



local energy emission, e m = L/m, Cp is the specific heat 
capacity at constant pressure, V is the radiative temperature 
gradient, V a( j is the adiabatic temperature gradient and V M 
is the mean molecular weight gr adient. For simp li city, the 
thermodynamic ratio j used by iMaeder fc Zahnl (| 19981 ) is 
set to unity. This is the correct for a perfect gas. We have 
also approximated the factor g/g of I Zahnl JT992) by 4 ^,^ . 
This is a suitable approximation throughout most of the 
star. The constant Co is included for aid of calibration and 
is discussed further in section [2/T1 

The same formulation for the meridional circulation is 
used i n all the test models except those that use the formal- 
ism of iHeger. Langer fc Wooslevl (2000) where the circula- 
tion is treated as a diffusive rather than an advective process 
and the characteristic circulation velocity is calculated from 
Eddington-Sweet circulation. 



2.4 Mass loss with rotation 

Observational evidence for enhanced mass loss is mixed (see 
IVardvalll985l ; lNieuwenhuiizen fc de JagedH9"88h but theoret- 
ically near-critical rotation must drive additional mass loss 
to remove angular momentum and prevent the surface of th e 
star from rotating super-critically (Friend & Abbott 1986). 
We use the enhanced mass-loss rate of lLangerl (1 19981 ) 



M = M n =o 



1 



(11) 



where we take £ = 0.45 and f2 crit = \J^§^r- 
A more complete discus sion of the critical rot ation 
rates of stars is given by iMaeder fc Mevnetl (|2000l ) and 
iGeorgy. Mevnet fc Maederl 1 (20 111 ). 



2.5 Rotation in convective zones 

Current ID models of stellar evolution generally assume that 
convective zones are kept in solid body rotation. This may be 
caused by str ong magnetic fields induced by dynamo action 
(|Spruidll999T ) but there is no conclusive evidence that real 
fields generated by this mechanism are strong enough to en- 
force solid body rotation. Certainly in the Sun we see latitu- 
dinal variations i n the angular velo city throughout the outer 
convective layer l|Schou et al.lll998l ). Standard mixing length 
theory suggests that a rising fluid parcel should conserve its 
angular momentum before mixing it with the surrounding 
material after rising a certain distance. This would lead to 
uniform specific angular momentum rather than solid body 
rotation. This is suppor ted by a recent MLT -based closure 
model for rotating stars jLesaffre et alj|20ich and by 3D hy- 
drodynamic simulations l|Arnett fc Mae kin 2009). In reality 
magnetic fields are likely to play some part in the transport 
of angular momentum but it is uncertain whether these are 
strong enough to affect the hydrodynamics. The asymptotic 
behaviour of the rotation profile in convective zones could 
have a profound effect on the evolution of the star, first 
because the total angular momentum content of a star for 
a given surface rotation increases dramatically for uniform 
specific angular momentum and secondly because uniform 
angular momentum in the convective zone produces a layer 
of strong shear at the boundary with the radiative zone and 
this drives additional chemical mixing. To explore the differ- 
ent possible behaviours we have introduced, in RoSE, the 
capacity to vary the distribution of angular momentum in 
convective zones as discussed in section [2~6l 



2.6 Angular momentum transport 

Differential rotation is expected to arise in stars because of 
hydrostatic structural evolution, mass loss and meridional 
circulation. Because of this stars are subject to a number 
of local hydrodynamic instabilities. These are expected to 
cause diffusion of radial and latitudinal variations in the an- 
gular velocity in order to bring the star back to solid body 
rotation, its lowest energy state. This occurs with char acter- 
istic d iffusion coefficients D s h ca r and respectively. IZahnl 
( 1992) proposed that, because of the strong stratification 
present in massive stars, the turbulent mixing caused by 
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these instabilities is much stronger horizontally than verti- 
cally (Dh S> -Dshear)- This leads to a situation where the 
angular velocity variations along isobars are negligible com- 
pared to vertical variations. Furthermore, all other state 
variables are assumed to be roughly constant over isobars 
and the mixing produces horizontal chemical homogeneity. 
This is referred to as shellular rotation, where we describe 
the angular velocity by Q. = f2(r). 

Taking into account all of the processes described in sec- 
tion^ we us e the evolution equation for the angular velocity 
l|Zahnll 19921 ) 



and 



d{r 2 fl) 

dt 



1 d(pr 4 nU) 1 d ( 4 dO. 
bpr 2 or pr £ or \ or 



+- 



1 d 



pr 2 dr 

and for the chemical evolution 

Oct _ 

dt r 2 dr 



pD con 



„(2+«) 



dr { 



2-n) 



(-Dshear + Aff + Ar2=o) T' 



Or 



.dci 



dr 



(12) 



(13) 



where Ci is the abundance of species i. The diffusion co- 
efficient -Dconv is non-zero only in convective zones and the 
D shear and -Deft are non-zero only in radiative zones. The pa- 
rameter n sets the steady state specific angular momentum 
distribution in convective zones, n — 2 corresponds to solid 
body rotation and n = corresponds to uniform specific 
angular momentum. The coefficient Aff describes the effec- 
tive diffusion of chemical elements because of the interaction 
between horizontal diffusion and meridional circulation so 

D« = &. (14) 



30 Ar 



T he main differences between models llTalon et al 



19971 ; iHeger. Langer fe Wooslevl |2000| ; iMevnet fc Maeder 
200(3 : iMaederl 120031 . e.g.) are the treatment of the merid- 



ional circulation, U, the diffusion coefficients, Ahcar, Aff 
and -Dconv, and the steady power-law distribution of angu- 
lar momentum in convective zones determined by n. We 
describe many of these models which we have implemented 
with RoSE in section |2~71 

2.7 Test cases 

In this paper we consider a number of common models for 
comparison. The details of the models are summarized in 
Table [1] Unless otherwise stated, the metallicity is taken to 
be solar (Z = 0.02). For each model we compute the stellar 
evolution for a range of masses between 3 M0 and 100 Mq 
and initial equatorial surface rotation velocities between 
and 600 kms -1 , except where the initial surface rotation rate 
would be super-critical. 

2.7.1 Case 1: 

Here we use the formulation for -Dshear o f lTalon et al.1l|l997h . 



A 



N 2 /{K + D^) + N 2 /D^ 



(15) 



where Ri c = (0.8836) 2 /2 is the critical Ric hardson numb er 
which we have taken to be the same as did lMaederl l|2003h . 

^ = -"tu(M p J^- v> (I6) 



fleff / <91n_p\ dlnp 
Hp \ d In fi J pT d In P 



Following [Maederi 120031 ) we take 

D h = 0.134r (rW [2V - aU])i 



where 



1 d(r 2 Q.) 



(17) 



(18) 



(19) 



The differential equations derived bv lZahnl 1^)92) are fourth 
order in space. Our model differs in that third order deriva- 
tives and above cannot be reliably computed and must be 
ignored. The constant Co is included as a means of calibrat- 
ing the model in light of this difference. Because our ultimate 
intention is to compare these models to data from the VLT- 
FLAMES survey of massive stars, we have chosen Co so 
that we reproduce the terminal-age main-sequence (TAMS) 
nitrogen enrichment of a 40 Mq star initially rotating at 
270 km s" 1 with galactic composition given by iBrott et al.l 
<|201ll ). This gives Co = 0.003. Whilst this is admittedly 
quite small, it is important to note that the non-linearity 
of the angular momentum transport equation means that a 
small change in the amount of diffusion corresponds to quite 
a large change in Co. We could have similarly adjusted the 
magnitude of D s hcar in case 1 to give the desired degree 
of nitrogen enrichment. However, attempting to reduce the 
diffusion coefficient results in more shear which drives the 
diffusion back up so Dshear would need to be reduced by a 
large factor to give the desired effect. This would have the 
further consequence that the meridional circulation would 
have a much stronger effect on the system and produce stars 
with a very high degree of differential rotation between the 
core and the envelope. Hence we have chosen to adopt this 
method of calibration as the most physically reasonable. 
The diffusion of angular momentum in convective zones is 
determined by the characteristic eddy viscosity from mixing 



length theory such that D c . 
model takes n — 2. 



D 



mlt — ^fmlt^mlt- 



This 



2.7.2 Case 2: 

This is the model of IHeger. Langer fc Wooslevl (|2000l ). In 
this case we set (7 = because circulation is treated as a 
diffusive process. The details of the various diffusion coef- 
ficients are extensive so we refer the reader to the original 
paper. With their notation the diffusion coefficients are 



Ashe 

and 



D>sem + D>dsi + Asm + D>ssi + Ass + D>gsf (20) 



Deff = (fc - 1)(D>dsi + Asm + Assi + Ass + -Dgsf), (21) 

where each of the A corresponds to a different hy- 
drodynamical instability. We note that our notation dif- 
fers slightly from the original paper. The diffusion co- 
efficie nts v and D used by IHeger. Langer fc Wooslevl 
(2000) are equivalent to Dshear a nd A hear + Ah respec- 
tively. IHeger. Langer fc Wooslevl l|2000D take f c = 1/30 
which is what we use here. T he parameter used in 
IHeger. Langer fc Wooslevl (|2000h is taken to be zero. Mean 
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Table 1. The diffusion coefficients used by each of the different models. 



Case 


n 


^shcar 


D h 


^conv 


1 


2 


Talon et al. C f 997) 


Maeder (2003) 


fmlt 


2 


2 


Heeer. Laneer & Wooslev (2000) 


N/A 


D m lt 


3 


2 


Zahn (19921 


Zahn (1992) 


D mlt 


4 





Talon et al. ( 1997) 


Maeder (2003) 


D m lt 


5 





Talon et al. (1997) 


Maeder (2003) 


10 10 cm a s _1 


6 





Talon et al. ( 1997) 


Maeder (2003) 


10 14 cm 2 s- 1 







molecular weight gradients play an important part in chemi- 
cal mixing near the core however we have performed a num- 
ber of test runs with / M = 0.05. Although there were some 
differences, they were not significant and could be largely 
masked by modifying the other free parameters associated 
with t his case. The model differs f rom t he original formula- 
tion of iHeger. Langer fc Wooslevl (|2000l ) in that we are un- 
able to use STARS to consistently compute non-local quan- 
tities such as the spatial extent of instability regions used in 
some of the expressions for the various diffusion coefficients. 
We do not expect this to have a significant effect on the re- 
sults since Des dominates the total diffusion coefficient and 
its limiting length scale is the pressure scale height rather 
than the extent of the unstable region. 

This model is calibrated by modifying the dominant 
diffusion coefficient for transport owing to meridional circu- 
lation, .Des, by a constant of order unity to give the same 
TAMS nitrogen enrichment as case 1 for a 20 Mq star with 
initial surface angular velocity of 300 km s" 1 . This model has 
n = 2 and D CO nv = Anit- 

2.7.3 Case 3: 

This is a reproduction of the original model of lZahnl (|l992T I 
and is included as a baseline to highlight the differences in 
predictions of stellar rotation from the original model. In 
this model 

-Dshcar = ^v,h + ^v,v, (22) 

where 

8C1 Ric v ( r dn Y 

K is the thermal diffusivity and a is the same as in case 1. 
In this model 

D h = C 2 r\2V - aU\, (25) 

and Ci and Ci are constants used for calibration. We con- 
strain C\ and C2 by matching as closely as possible the 
TAMS nitrogen enrichment and luminosity of a 20 Mq case- 
1 star with initial surface rotation of 300kms _1 . We find 
that Ci = 0.019 which is surprisingly small. This is because 
this model does not take into account mean molecular weight 
gradients and this leads to far more mixing between the con- 
vective core and the radiative envelope than in case 1. The 
TAMS luminosity is always greater than case 1 and so we 
minimize it with respect to C2 so C2 — ► 00. This is realised 



by setting D a a = and is the case when the horizontal dif- 
fusion completely dominates over the meridional circulation 
and so is consistent with our assumption of shellular rota- 
tion. 

Apart from ignoring mean molecular weight gradients, 
the main objection to this model is that in the formulation of 
Dh we assume that, if the horizontal variation in the angu- 
lar velocity along isobars takes the form Cl = S72(r)P2(cos6 l ), 
then Q%(r) /Cl(r) is constant and this is not physically justi- 
fied. Again we set n = 2 and D conv = D m i t . 

2.7.4 Cases 4, 5 and 6: 

For these cases we use the same -Dahear and Dh as in case 1 
but we set n = to produce uniform specific angular mo- 
mentum through the convective zones and test the effect of 
varying the convective diffusion coefficient D conv . 

fAnit case 4, 

-D conv = < 10 10 cm 2 s" 1 case 5, (26) 
[ 10 14 cm 2 s -1 case 6. 

2.7.5 Metallicity variation 

We have also calculated the evolution in cases 1, 2 and 3 
for the same masses and velocities but with Z = 0.001. We 
shall represent these cases with a superscript Z (e.g. Case 
l z is the low metallicity analogue of Case 1). 



3 RESULTS 

Whilst there are many potential observables which may be 
used to distinguish between different models, it is impor- 
tant to choose the ones that are most easily compared with 
observational data. From our stellar evolution calculations 
we find a number of important differences between the test 
cases. 

3.1 Evolution of a 20 Mq star in cases 1, 2 and 3 

First it is helpful to briefly examine the internal evolution 
that occurs. We consider here the main-sequence evolution 
of a 20 Mq star with initial surface rotation of 300 kms -1 for 
cases 1, 2 and 3. Although the centrifugal force causes some 
change in a star's structure, its evolution is strongly affected 
by changes in the chemical composition. Fig[T]shows how the 
composition of the rotating 20 Mq case-1 star differs from 
a non-rotating 20 Mq star at the end of the main sequence. 
The difference in the rotation-induced mixing produces the 
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Figure 1. Terminal-age main sequence composition of 20 Mq case-1 stars. The solid lines are for a star initially rotating at 300 kms 
and the dashed lines are for a non-rotating star. Both stars have central hydrogen abundance Xh = 10~ 3 . Note that rotational mixing 
results in a larger core and mixing of helium and nitrogen throughout the radiative envelope. 



variation in results between the various cases. Fig. [2] shows 
the angular velocity profile and the diffusion coefficient for 
vertical angular momentum transport in radiative zones pre- 
dicted by each of cases 1, 2 and 3 at the zero-age main se- 
quence (ZAMS). Note that, even though the stars have the 
same surface rotation, their core rotation and hence total 
angular momentum content can vary significantly between 
models 

Despite their similar treatments, cases 1 and 3 have 
quite different initial rotation profiles. This is largely due to 
our choice of calibration. Because in case 3 we ignore mean 
molecular weight gradients, the overall efficiency of mixing 
must be reduced to match our calibration criterion (section 
12.70 . This means that shear diffusion is much weaker relative 
to advection and so a profile with more differential rotation 
results. Had we chosen to calibrate the mixing by reducing 
Co instead of Ci we would have found the opposite effect. 
This highlights one possible pitfall of including multiple free 
parameters within a given system. 

Case 2 is dominated by diffusion because of the diffusive 
treatment of meridional circulation. Recall that meridional 
circulation is treated advectively in cases 1 and 3 so is not 
included in the diffusion coefficient. In fact it is responsible 
for production of the shear at the ZAMS despite turbulence 
trying to restore solid body rotation. Because there is no 
perturbation to the rotation at the start of the main se- 
quence, the star in case 2 rotates as a solid body. As the 
star evolves and mass is lost from its surface the solid body 
rotation is disturbed. Even so, because of the strong diffu- 



sion, case 2 stars never deviate far from solid body rotation 
as can be seen in Fig. [3] We note that case-1 stars reach the 
end of the main sequence with a higher mass than those in 
case 2 or case 3. This is because case-2 and case-3 stars have 
a longer main-sequence life owing to more efficient mixing 
at the core-envelope boundary. This allows hydrogen to be 
mixed into the core more rapidly than in case 1. This also 
leads to larger core masses in cases 2 and 3 compared to 
case 1. 

We see in Fig.[2]that the predicted diffusion coefficients 
for cases 1 and 2 are similar throughout most of the en- 
velope. By the TAMS, the diffusion predicted by case 2 is 
significantly lower than the other two cases. This is possibly 
because rising shear owing to rapid hydrostatic evolution at 
the TAMS causes the diffusion in cases 1 and 3 to increase 
while in case 2 diffusion is dominated instead by the circu- 
lation. Unsurprisingly, the diffusion coefficient in case 3 is 
similar in form to case 1 but significantly smaller, a result of 
our choice of Ci. We note however that the diffusion coeffi- 
cient predicted by the two cases is very similar by the end 
of the main sequence. Also, the diffusion coefficient at the 
core-envelope boundary at the end of the main sequence in 
case 1 is around an order of magnitude lower than in case 3. 
This is partially due to the core of the case-3 star rotat- 
ing faster than in case 1 but is mostly due to the inclusion 
of the me an molecular weigh t grad ient in the formulation 
for case l. lFrischknecht et alJ l|20ld ) (who use a model very 
similar to our case 1) predict a much greater decline in the 
mixing near the core but we have been unable to repro- 
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Figure 2. Zero-age main sequence properties of a 20 Mq star initially rotating at 300 km s 1 . The left panel shows the angular velocity 
through the star and the right panel shows the diffusion coefficient for chemical transport through the radiative envelope 
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Figure 3. Terminal-age main sequence properties of a 20 Mq star initially rotating at 300kms — 1 . The left panel shows the angular 
velocity through the star and the right panel shows the diffusion coefficient for chemical transport through the radiative envelope. The 
vertical line of the left panel separates the convective and radiative zones of the case-1 star. 



duce this. It is likely that the difference is a result of the 
inhibiting effect of the mean molecular weight gradient on 
the rotation al mixing. Whilst it is inc luded in our models, 
the results of lFrischknecht et al.l |2010t ) suggest that towards 
the end of the main-sequence its effect covers a much larger 
proportion of the radiative envelope than in our models. 



3.2 Effect on HR diagram 

As expected, the effects of rotation on the structure and 
chemical evolution of each star are significant across the HR 
diagram. We have plotted the case 1 models in Fig [4] Cen- 
trifugal forces cause the star to expand making it dimmer 
and redder. However, because of the additional chemical 
mixing, more hydrogen is made available during the main 
sequence so, by the main-sequence turn off, rotating stars 
are generally more luminous than similar mass non-rotating 
stars. This effect becomes more pronounced at higher rota- 
tion rates and is most apparent for stars with masses up to 
20 M@. 



There are clear differences between the predicted evolu- 
tion of rotating stars in each of the three cases. Fig.[5]shows 
the evolution of five different stellar masses in the HR di- 
agram for cases 1, 2 and 3 with an initial surface rotation 
of 300 km s" 1 . These cases are the three different models of 
rotational mixing in radiative zones at solar metallicity. Ro- 
tating case-1 stars appear to be the most luminous at low 
masses but least luminous at high masses. Case-2 stars are 
also consistently cooler than their case-1 and case-3 equiv- 
alents except below 10 Mq. This is because the strength of 
rotation-induced mixing increases rapidly with mass in case- 
2 stars unlike in cases 1 and 3 where the difference is more 
modest. 

Although there is apparently a large difference between 
the three models in the HR diagram, to distinguish be- 
tween them from stellar populations requires either a very 
large sample or accurate independent measurements of stel- 
lar masses and rotation velocities. Both of these are very 
challenging but, with the advent of large scale surveys, the 
former is quickly becoming a possibility. 
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Figure 4. Stellar evolution tracks for stars between 5 Mq and 100 Mq calculated with RoSE for case 1 with surface velocities of 
to 400 km s" 1 . The tracks for cases 4, 5 and 6, where the angular momentum transport in convective zones is varied, are almost 
indistinguishable from each other but produce more luminous stars than in case 1. There is significant difference between the predictions 
made for cases 2 and 3, where the angular momentum transport in radiative zones is varied. 



3.3 Nitrogen enrichment 

Currently the key test for any model of stellar rotation 
is how w ell it reproduces th e spread of data in a Hunter 
diagram (|Hunter et al.ll2009l ). Hunter plotted nitrogen en- 
richment against surface rotation. Large scale surveys, suc h 
as the VLT-FLAMES Tarantula survey l|Evans et ai.ll2010l ). 
will greatly increase the data available for surface rotation 
velocities and surface abundances over the coming decade. 
Thus, if different models can be distinguished in a Hunter 
diagram, this would form a key test for stellar rotation mod- 
els. 

In Fig. [6] we plot our theoretical Hunter diagram for 
10 Mq and 60 Mq stars at different initial surface rotation 
velocities with the different radiative-zone models at solar 
metallicity. Each star begins at the bottom of the plot with 
the same nitrogen abundance but different initial surface 
velocities. There is an initial period where the star spins 
down before any enrichment has occurred. During this time 
the star moves straight to the left of the plot. Eventually 
the surface nitrogen abundance begins to increase because 
of rotation induced mixing from the burning region. At the 
same time the star continues to spin down because of mass 
loss and structural evolution. The net effect is that the star 
moves towards the upper left-hand region of the plot. At 
the end of the main sequence the star expands and rapidly 
spins down. This appears as a near-horizontal line as the 



star moves rapidly towards the left-hand edge of the plot. 
Some further enrichment may occur during the giant phase. 

The evolution of the surface abundances is very model 
dependent. Case 3, which is based on the early model of 
IZahnl (|l992T l. gives more nitrogen enrichment than case 1 
at high masses (M ^ 60 Mq) but significantly less at low 
masses (M ?C 10 Mq). Most notably though, the case-3 stars 
spin down to a far greater degree before enrichment occurs. 
We attribute this to the neglect of mean molecular weight 
gradients. Mixing near the core is more efficient in case 3 but, 
due to the overall calibration, is weaker near the surface. 

For case 2 the amount of mixing in lower-mass stars 
is much less than for both cases 1 and 3. By comparison 
the enrichment of case-2 60 Mq stars is greater than in the 
other two cases, particularly for slower rotators. This mass 
dependent behaviour of the rotating models could provide 
important clues to distinguish between the models. 

At solar metallicity, owing to the enhanced mass loss, 
massive stars spin down before the end of the main sequence 
so, in this case, we would not observe the absence of the mod- 
erat ely rotating, highly enriched stars seen at low metallic- 
ity (|Hunter et alj I2009T ) . Observations of multiple clusters 
at different ages at this metallicity would be a good test 
for rotating stellar models because the evolution across the 
Hunter diagram is significantly different in each case even 
when they are calibrated to give the same level of enrich- 
ment at the end of the main sequence. 
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Figure 5. St ellar evolution tracks for s t ars b etween 5Mq and 100 Mq with initial surface rotation of 300 km s -1 for cases 1, 2 and 3. 
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Figure 6. Nitrogen enrichment (by number of nuclei) variation with initial surface rotation for cases 1, 2 and 3. Stars start on the ZAMS 
with low nitrogen abundances and high velocities and evolve to higher abundances and lower velocities during the main sequence. The 
left panel is for 10 Mq stars and the right panel is for 60 Mq stars. Note the different scales for each panel. As expected the enrichment 
is much greater for more massive stars. Case-3 stars spin down more before any enrichment occurs. They give greater enrichment than 
either of the other two models at high masses but significantly less than case 1 for low masses. Cases 1 and 2 enrich to a similar degree 
for the high-mass stars but case 2 exhibits significantly less enrichment for low-mass stars than case 1. 



3.4 Helium-3 enrichment 

Apart from the enrichment of nitrogen, rotation can have a 
large effect on the evolution of other elements. Changes in 
the carbon and oxygen abundances in rotating stars pre- 
dicted from models have been considered but the accu- 



racy of th e data prohibits any stro ng conclusions from be- 
ing made. iFrischknecht et all (|2010t ) discuss the effect rota- 
tion may have on the surface abundances of light elements. 
We consider here the evolution of the surface abundance 
of He. A similar analysis could be performed for other 
elements such as lithium and boron. The changes in the 



© 0000 RAS, MNRAS 000, 000-000 



10 A. T. Potter, C. A. Tout & J. J. Eldridge 



overall abundance of 3 He because of rotation could partially 
explain the discrepancy between the predicted abundances 
produced by stars and the lack of enrichment of the ISM 
co mpared to levels predicted by primordial nucleosynthe- 
sis dDearborn. Schramm fc Steigm an 1986; Hata et all 19951 ; 
iDearborn. Steigman fc Tosilll996h . Th is has been exp lained 
in the past by thermohaline mixing (|Stancliffdl2O10h but, 
as our results show, the surface 3 He abundance is strongly 
affected by rotation and so it is likely to make at least some 
contribution to this effect. We leave the issue of whether the 
total production increases or decreases over the stellar life- 
time to future work. In either case, the evolution of helium-3 
with respect to the surface rotation is very different between 
alternative models and so, as for nitrogen enrichment, this 
would form a useful means to compare stellar rotation mod- 
els. Unlike nitrogen, helium-3 enrichment is stronger at low 
masses and so provides a greater number of candidate stars 
for observations. 

Fig. [7] shows the helium-3 enrichment for 10 Mq and 
60 Mq stars of varying initial surface velocities for each of 
the different radiative zone models at solar metallicity. As 
for nitrogen, all three cases show comparable amounts of 
enrichment in high-mass stars. The amount of enrichment 
at the end of the main sequence is the same in each case but 
case-3 stars are slightly more enriched at all rotation rates. 
Case-3 stars spin down much more before enrichment occurs 
so the paths for these stars lie to the left of case 1 and 
case 2 stars but the amount of enrichment at the end of the 
main sequence is comparable, though slightly higher than, 
the other two cases. 

The difference between the test cases is far greater at 
lower masses. Both case-2 and case-3 stars show substan- 
tially less enrichment during the main sequence especially 
for slow rotators and case-3 stars have much slower surface 
rotation at the end of the main sequence than in the other 
two cases. Both of these contribute to very different evolu- 
tion which should be distinguishable observationally. Indeed, 
whilst the models may produce similar results for full pop- 
ulation synthesis calculations, it has been found that there 
is often far less agree ment when d ifferent mass ranges are 
considered separately (|Brott et al. II2OI0T ). 

3.5 Metallicity dependence 

In order to compare stellar rotation models with data it 
is important to distinguish which effects are observable at 
different metallicities. Low-metallicity stars are particularly 
useful because they have s i gni ficantly lower mass-loss rates 
ijVink. de Koter fc Lamersl l200lf ). For stars of metallicity 
Z = 0.001 the mass-loss rate is roughly ten times lower 
than at solar metallicity. This allows us to rule out effects 
on the models produced by our prescription for mass loss. 

The low-metallicity cases show similar distinctions in 
the HR diagram to those at solar metallicity, although 
high-mass, rapidly rotating, case-2 stars are sufficiently well 
mixed to undergo quasi-homogeneous evolution. There are 
also significant differences in the nitrogen enrichment be- 
tween the different models. Because the mass-loss rate is 
lower in low-metallicity stars they retain their surface ro- 
tation for much longer so the main sequence appears much 
more vertical in a Hunter diagram. Unlike the solar metal- 
licity cases, case 2 Z exhibits significantly more mixing than 



case l z particularly for slow and moderately rotating stars 
(Fig. |8]) . This is the complete opposite of the results at solar 
metallicity and highlights the importance of testing different 
stellar environments to discover clues to distinguish between 
different models. In contrast, the enrichment of helium-3 in 
case 2 Z stars is less than in case l z . 



3.6 Surface gravity cut-off 

As a consequence of increasing stellar radius and angular 
momentum conservation all models for stellar rotation pre- 
dict a rapid decay in the surface rotation velocity after 
the end of the main sequence. Observations suggest that, 
even for rotating stars, there is a sharp cut-off in the effec- 
tive gravity at log 10 (g/cm 2 s -1 ) ~ 3.2 when a star leaves 
the main sequen ce and moves over to the giant branch 
(IBrott et al . 2011). This effect depends on stars reaching the 
TAMS without spinning down too rapidly during the main 
sequence. Therefore it is more easily seen at lower metallic- 
ities where the mass-loss rate is lower. The observed value 
for the TAMS gravity can be enforced in rotating models 
by including a degree of overshooting. However this simply 
introduces an additional free parameter into the models. In 
Fig. [5] we show the different cut-offs in the effective grav- 
ity predicted by cases l z , 2 Z and 3 Z . The end of the main 
sequence is indicated by a distinct cusp in the path of the 
star in the range 3 < log 10 (<;/cm 2 s _1 ) < 4. Note that the 
expected number of stars with significant rotation after the 
main-sequence cut-off is low because the evolution to a very 
slowly rotating giant is extremely rapid compared to the rest 
of the main sequence. 

For higher-mass stars there is a sharp cut-off in the sur- 
face gravity at the end of the main sequence that is the same 
regardless of the model used although the TAMS surface 
rotation is somewhat higher in case 2. The mixing in low- 
metallicity, high-mass, rapidly rotating, case-2 stars is very 
efficient and so they evolve almost homogeneously and thus 
they appear differently in the plot. For lower-mass stars the 
change in surface gravity at the end of the main sequence 
is still clear but varies by around an order of magnitude 
across all rotation rates and test cases. Case-2 and case-3 
stars have lower terminal-age main-sequence surface grav- 
ities than case-1 stars and show a tendency towards lower 
TAMS surface gravities for more rapid rotators. Case-3 stars 
generally have lower rotation rates at the end of the main 
sequence for low-mass stars. This distinguishes them from 
case 2. Once again we see that the difference in stellar prop- 
erties predicted by each model is strongly dependent on rota- 
tion rate, mass and metallicity suggesting that it is essential 
to explore populations covering as wide a range as possible 
in order to test rotating models. 

3.7 Alternative models for convection 

In order to compare the difference in the evolution of a star 
owing to the details of the model for convective angular mo- 
mentum transport we now focus on cases 1, 4, 5 and 6. Uni- 
form specific angular momentum in the core causes more 
shear mixing near the core-envelope boundary than when 
the core is solid body rotating as shown in Fig. [10] This re- 
sults in higher luminosity stars with similar temperatures. 
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Figure 7. Helium-3 enrichment variation with initial surface rotation for cases 1, 2 and 3. The left panel is for 10 Mq stars and the 
right panel is for 60 Mq stars. At high masses all three cases show a similar degree of enrichment though there is some variation in the 
amount. Case-3 stars spin down more before enrichment occurs and so lie to the left of the other two cases. At low masses, both case-2 
and case-3 stars are less enriched at the TAMS than case-1 stars. This is especially true for slow rotators. 




Figure 8. Chemical enrichment variation at low metallicity. The left panel shows the enrichment of nitrogen and the right panel shows 
the enrichment of helium-3. The left-hand plot is for 60 Mq stars and the right-hand plot is for 10 Mq stars. Case 2 Z now shows much 
higher nitrogen enrichment than case l z , the opposite to what we found in the solar metallicity cases. The enrichment of helium-3 in 
case 1 Z stars is still considerably less than in case \ z . 



There is almost no difference between cases 4, 5 and 6 in the 
HR diagram. 

When we compare the different models for convection 
in a Hunter diagram we see that cases 4, 5 and 6, which 
have uniform specific angular momentum throughout their 
convective zones, have significantly more enrichment for all 
masses and rotation rates than case 1. The difference is more 
pronounced for higher-mass rapid rotators (Fig. How- 
ever, we note that it is more difficult to distinguish between 
cases 4, 5 and 6. For the highest mass stars we do find some 
difference in the enrichment of nitrogen and helium-3 be- 
tween the models but recall that there is a difference in 
fl conv of four orders of magnitude between cases 5 and 6. 

Given the small magnitude of the change in enrichment 
and structure over such a range of diffusion coefficients it 
seems unlikely that these tests can adequately distinguish 
between convective models. In addition, adjusting the cali- 
bration of case 1 could produce a very similar effect making 
it difficult even to distinguish between n = and n — 2 



models from observations. At the same time though, it is 
interesting to note the significant change that modifying the 
core angular momentum distribution has had on the evolu- 
tion of the surface composition. 



4 CONCLUSIONS 

Rotation in stars has a number of profound effects on their 
evolution. Not only are there significant changes in the hy- 
drostatic structure |Endal fc Sofialll976l ) but this causes a 
thermal imbalanc e that can le ad to a strong meridional cir- 
culation current (ISwee tl ll950h . Meridional circulation leads 
to additional shear which induces a number of instabilities. 
The resulting turbulence leads to strong mixing of both an- 
gular momentum and chemical composition. 

Although most modellers include all of these effects, 
the exact implementa tion of stellar rotation can vary dra- 
matically. For example lHeger. Langer fc Wooslevl |2000) use 
a model where meridional circulation is treated diffusively. 
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Figure 9. The effective gravity variation with surface velocity for 10 Mq stars (left panel) and 60 Mq stars (right panel) for cases l z , 
2 Z and 3 Z . The end of the main sequence occurs at the right-most cusp of each path. Although the surface gravity is similar for all the 
models and all initial rotation rates, the TAMS surface velocities are very different between the different cases. Exceptions to this are 
the rapidly-rotating, high-mass, case-2 z stars which undergo almost homogeneous evolution. 




Figure 10. Angular velocity distributions for 60 Mq stars initially rotating at 300 km s^ 1 for different convective models. The left panel 
shows the ZAMS distributions and the right panel shows the TAMS distributions. Despite there being four orders of magnitude difference 
between the convective diffusion coefficients in cases 4,5 and 6, there is very little difference between the angular velocity distributions 
they predict. In each case, the models with uniform specific angular momentum in convective zones predict more shear at the convective 
boundary than the models with solid body rotation in convective zones. 



The diffusion owing to shear is the linear combination of a 
number of coefficients based on different poss ible instabil i- 
ties. On the other h and, models such as those o f lZahnl (19921 ). 
iTalon et al.1 j 19971 ) and lMaederl \200$ ) treat the circulation 
as advective and use a single diffusion coefficient based on 
the magnitude of Kelvin-Helmholtz instabilities induced by 
shear. In these models it is also necessary to define the mag- 
nitude of diffusion along isobars and this too varies between 
different models. 

We have shown that different models generally give rise 
to similar qualitative conclusions but there are significant 
differences in the results based on mass, rotation rate and 
metallicity. There are also open questions about how angular 
momentum transport occurs in convective zones. 

Compar ing th e mo dels based on ITalon et al.l 

il 19971) and iMaederl (l2003h . cas e 1, and that based on 
iHeger. Langer fc Woosle^l^OOOT) , case 2, we find that 
case 1 gives higher luminosity stars for masses less massive 



than 10 Mq and more luminous, hotter stars at higher 
masses than case 2. High-mass stars give similar levels of 
nitrogen enrichment in each case but case 2 produces far 
less enrichment for low-mass and intermediate-mass stars 
than case 1. The situation is similar for their helium-3 
enrichment. 

The predicted effects of rotation appear to be highly 
dependent on the metallicity. This is one of the clearest 
tests for different stellar rotation models. At metallicity of 
Z — 0.001, case 2 actually produces significantly more ni- 
trogen enrichment in high-mass stars although the degree of 
helium-3 enrichment is still lower in case 2 than case 1. Ad- 
ditional effects are seen in lower-metallicity stars where the 
mass-loss rate is lower such as the variation of the surface 
gravity with respect to surface rotation velocity. We see a 
sharp cut-off in the effective gravity at the end of the main 
sequence but the TAMS rotation rates are very different 
between different models. Case-2 stars reach the end of the 
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Figure 11. Nitrogen enrichment variation with initial surface rotation for cases 4, 5 and 6. The left panel is for 10 Mq stars and the right 
panel is for 60 Mq stars. As expected the enrichment is much greater for more massive stars. There is much more enrichment for models 
in which the core has uniform specific angular momentum, an effect that becomes progressively more pronounced for higher masses and 
rotation rates 



main sequence with higher rotation rates than the other two 
cases for both low and high-mass stars. Case-1 stars reach 
the end of the main sequence with higher surface gravity 
than the other two cases but only for lower-mass stars. 

All current models treat convective zones as rotating 
solid bodies. This may be justified if convecti ve zones can 
generate sufficiently strong magnetic fields fe.g. lSpruitlll999l ) 
but, if not, hydrodynamic models and calculations suggest 
that convective zone s should tend towards uniform specific 
angu lar momentum (|Lesaffre et alj |2010| ; lArnett fc Maekinl 
l2009h . Identifying whether this is the case or not is difficult 
from surface observations. There is no significant change in 
the paths of stars across the HR diagram between cases 4, 5 
and 6. In fact, there is no test we have found that would ade- 
quately differentiate between these three cases even though 
the diffusion coefficient in convective zones varies by four 
orders of magnitude. There are some minor differences in 
the amount of enrichment for massive stars but these would 
be hard to test with existing data. The models with uni- 
form specific angular momentum generally produce slightly 
more luminous stars and higher surface chemical enrichment 
than those models with solid body rotating cores. However, 
slightly less efficient mixing in the radiative zone could mask 
this difference. 

We have thus far not included magnetic fields in the 
models. It has been suggested that the strong turbulence 
generate d by rotation could result in a radiative magnetic 
dynamo (|Spruitl 119991 ). A sufficiently strong magnetic field 
can effectively suppress the meridional circulatio n and re- 
duce the overall shear l|Maeder fc Mevnetl 12005). It could 
also result in additional mass and angular momentum loss 
(e.g. lLau. Potter fc Tout! libTll ) . As with rotation, there is 
little consensus on the details of magnetic field generation 
but it is generally accepted that the effects of rotation and 
magnetic fields cannot be considered in isolation. In the fu- 
ture we plan to include magnetic fields in RoSE in a similar 
manner in order to better explore this important feature of 
stellar evolution. 

Owing to the range of available models, it is an ex- 
tremely challenging problem to try to identify the one which 



best fits observed stellar populations. Now that large scale 
surveys are starting to produce data for many stars in dif- 
ferent regions it is becoming possible to make progress and 
isolate which effects dominate. The key to distinguishing the 
relevant physics seems to be taking measurements of groups 
of stars at different masses, ages and metallicities. In indi- 
vidual clusters, models should be able to match not only 
the full distribution of observed stars but also the expected 
distribution in each mass range. Whilst most of the models 
can be calibrated to fit the data for a single cluster, we have 
shown that the behaviour of each model is highly depen- 
dent on mass and metallicity and so being able to fit data 
for a range of stellar environments is the true test for any 
model. Using our new rotating stellar evolutio n code RoSE 
combined with populat ion synthesis codes (e.g. Ilzzard et al.l 
2009|; iBrott et al.ll201ll ) we plan on making a detailed com- 
parison of these models with the available data to identify 
how best to model stellar rotation. 
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